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Abstract. We find a close correspondence between certain partition functions of 
ideal quantum gases and certain symmetric polynomials. Due to this correspondence 
it can be shown that a number of thermodynamic identities which have recently 
been considered are essentially of combinatorical origin and known for a long time as 
theorems on symmetric polynomials. For example, a recurrence relation for partition 
functions appearing in the textbook of P. Landsberg is nothing else but Newton's 
identity in disguised form. Conversely, a certain theorem on symmetric polynomials 
translates into a new and unexpected relation between fermionic and bosonic partition 
functions, which can be used to express the former by means of the latter and vice 
versa. 
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1. Introduction 

The theory of ideal quantum gases, although a standard section in every modern 
textbook on statistical mechanics, is far from being a closed field of research. On one 
hand, there are continuous attempts towards a deeper understanding of the foundations 
of the field and to disclose unexpected regularities and relations between known facts. 
As examples we quote the work of M. H. Lee |^, on a unified treatment of Fermi and 
Bose gases and our earlier paper on a fermion-boson symmetry for harmonic oscillators of 
odd space dimensions 0]. Also this article is intended to contribute to such foundational 
questions 

On the other hand, there are various applications of the theory of ideal gases in 
almost all areas of modern physics, despite the apparent over-simplification of non- 
interacting particles. In nuclear physics the shell model is a prominent representative, 
which at low excitation energies describes the motion of noninteracting fermions in a 
common (mean) field 0. In solid state physics of one space dimension interacting 
fermions may be described as noninteracting bose gases (Luttinger liquids) [p, |^, |^. 
The crucial point is that the specific heat of fermions and bosons is the same for 
equidistant energy spectra, for recents works see p!0|] . 

In atomic physics the concept of ideal quantum gases helped to understand 



experiments on magnetically trapped atomic vapours 0, Here the system 

can be well described as an ideal quantum gas contained in an external harmonic 
oscillator potential, for an overview see ||T5| . Also the case of constant and small particle 
number N turns out to be of interest, thereby excluding the use of a grand canonical 
ensemble {N fluctuating) and/or the thermodynamic limit (A^ —>■ oo). Examples of 
applications, except the above-mentioned trapped systems, are quantum dots and 



fluorescence from a few electrons [l^] • For a comparison of canonical and grand canonical 



results in a two-state system see also 

In the context of these small quantum systems it is sometimes of interest to 
explicitely calculate the A^— particle partition function Z^, where -|- stands for bosons 
and — for fermions. For this purpose a certain recurrence relation for has recently 



been used by several authors [IS, 20, 0, H, 17]. It seems that none of the authors 



including us, was aware that this recurrence relation could be traced back to the time- 



honoured textbook of P. Landsberg |21|. Even more remarkable we find that, in some 
sense, Landsberg's formula was already discovered by Isaac Newton. Of course, Newton 
did not know anything about partition functions and Fermi/Bose statistics. However, 
there is an intimate connection between partition functions and symmetric polynomials, 
which allows to translate certain theorems on symmetric functions into statements about 
physical partition functions and vice versa. This will be the main subject of the present 
article. 

The basic idea will be spelled out in section 2. Since we write this article for people 
who may be experts in either field, combinatorics and statistical mechanics, but not 
in both, we will give some rather elementary introductions into the field of symmetric 
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polynomials (section 3) and partition functions (section 4). Section 4 also contains 
the "dictionary" needed for the translation between statements of the respective fields. 
Thus we obtain the following results, which now appear obvious, but have not yet been 
explicitely mentioned in the physics literature, namely that may be expressed by 
the Z~ ,n < N and vice versa. The most compact expression of this is the statement 
that the grand canonical partition functions satisfy Z~^{z)Z~{—z) = 1. In section 5 we 
add some closing remarks. 

2. The basic idea 

Before going into details we will describe the basic idea of the correspondence between 
partition functions and symmetric polynomials and illustrate this idea for the simple 
example of a system with L = 3 energy levels and N = 2 particles. 

It is well-known that the eigenstates of the A^— particle Hamiltonian (+ for 
bosons, — for fermions) of N non-interacting particles can be labelled by occupation 
number sequences {ni,n2, ■ ■ ■ jTIl), where JT-^ = N and n£ — 0,1 for fermions. 

The latter condition expresses the Pauli principle. Equivalently, the eigenstates are 
in 1 : 1-correspondence to the monomials x"^ x/-^'^ . . . x'^^^ of degree A^, where the 
Xi,X2, ■ ■ ■ ,Xl are (symbolic) commuting variables. In our example, wc obtain the mono- 
mials X1X2, xiXs, X2XS for the fermionic case, and, additionally, q, for the bosonic 
case. 

The eigenvalue of corresponding to the eigenstate labelled by 
is given by X^^^^ E^ne, if Ee denotes the i-th energy eigenvalue of the one-particle- 
Hamiltonian. Hence the corresponding eigenvalue of exp(— where, as usual, 
(3 — denotes the inverse temperature, will be 



This eigenvalue is nothing else but the value of the monomial ^ evaluated at 

Xi — e~^^^, . . . ,xl — e"^^^. To obtain the trace Trexp{—f3H^) we have to sum over all 
eigenstates, or, equivalently, over all monomials of degree N (subject to the constraint 
— 0,1 in the fermionic case). Interchanging summation and evaluation, we may first 
consider the bosonic polynomial (resp. the fermionic /jv) obtained by summing up 
all monomials of degree (satisfying the above constraint in the case of /jv) and then, 
in a second step, evaluate this polynomial at Xe = e~^^' ,1 = 1, . . . , L in order to obtain 
Tr exp(— /3i/^). The latter is usually called the partition function of the system. In our 
example, h2 = x\ + x\ + x\ + X1X2 + X1X3 + X2X-i and /2 = X1X2 + X1X3 + X2X3. Both 
polynomials are symmetric w. r. t. any permutation of their arguments xi,X2,X3. This 
also holds in the general case. 



L 




(1) 
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Thus we have estabUshed a connection between systems of N non- interacting bosons 
(resp. fermions) and certain symmetric polynomials 67V (resp. /jv). Specification of the 
one-particle Hamiltonian then corresponds to an evaluation of these polynomials at 
— e~^^^,£ = 1, . . . ,L and yields the particular partition function. 



3. Symmetrical polynomials 

We collect some well-known definitions and results on symmetric polynomials which are 
relevant for our purposes. Throughout we consider polynomials in a finite number of 
(symbolic) variables Xi, . . . Most results also hold for the case of infinitely many 
variables {n = 00) which is of physical interest, where the polynomials become formal 
power series. We will only indicate those cases where n = 00 becomes problematic. 
A polynomial p is called symmetric iff p{xi,...,Xn) = p{x„(^i), . . . ,Xa(n)) for all 
permutations u G The elementary symmetric polynomials fm can be defined 

through their generating function 

n n 

n(a: + x,) = ^/^x"-"^, m = 0,...,n. (2) 

e=i m=0 

For example, 

/o = 1, (3) 

/l ^ Xi + X2 + . . . + Xn, (4) 

/2 = ^XiXj, (5) 
i<j 

fn ^ X1X2 ■ ■ ■ Xn- (6) 

The letter "f ' for elementary symmetric polynomials is unusual, but chosen here to stress 
the association to "fermions". One important property of the /'s is that (fx, . . . , /„) 
forms a basis in the sense that any symmetric polynomial can be uniquely written as a 
polynomial of the For example, if n = 2, 

xl + xl^ {xi + X2f - 2xiX2 = /i - 2/2. (7) 

Of course, there exist other bases of symmetric polynomials. For our purposes, the 
following ones are the most important: 

Let be denote the sum of all monomials of degree £, for example 

bo = 1, (8) 

bi ^ Xi + X2 + . . . + Xn, (9) 

&2 = ^XiXj, (10) 
i<j 
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Obviously the 6^ are symmetric polynomials, sometimes called complete symmetric 
polynomials (our "b" refers to "bosons"). Also the sequence {bo, ... , hn) forms a basis in 
the ring of all symmetric polynomials in n variables, as we will see below. For example, 
if n = 2, 

xl + xl = 2{x\ + + X2X2) - {xi + X2)^ = 262 - bl- (11) 
Another basis is given by the power sums 

n 

s, = J2<' i = 0,...,n. (12) 

i=l 

(so = n has to be redefined in the case n = 00.) In order to prove that (sq, . . . , s^) 
forms a basis, it would suffice to show that every elementary symmetric polynomial fm 
could be expressed as a polynomial of the S£. This, in turn, follows immediately from 
the Newton identities , as they are called in Ref . P2| , 

»n— 1 

mU = Y,i-'^)'^'fkSm-k, l<m<n. (13) 

These identities could be viewed as recursion relations for the fm- Iterative substitution 
then yields the fm solely in terms of the sg, for example 

/3 = ^(S?-3S1S2 + 2S3). (14) 

Explicit formulae for these representations will be given below. 

Now consider the case of the bm- Here we have analogous identities which 
immediately imply that {bo, ... , 6„) will form a basis: 

m 

Y,^-lffkhm-k = 0, l<m<n. (15) 

fc=0 

These identities can most compactly be written in terms of the generating functions 

F{z)=Y,frZ^ = \{{l + X£z), (16) 

00 00 

i?(z) = 5^M^' = J](l-x,z)-i (17) 



r=0 



as 



B{z)F{-z) = 1. (18) 

Finally we give the explicit representations of the fm and bm in terms of the power 
sums se, as well as the representations of the fm in terms of the bk and vice versa. 
To write these representations in compact form, it is convenient to use partitions. A 
partition A of a positive integer A^, written as A h A^, is a way to write A^ as a sum, 
irrespective of order, in other words, a non-increasing sequence of positive integers 
(Ai, A2, . . . , Ap) such that Yli=i '^j — ^- = P is called the length of a partition. 

Another notation for A is (I™i2™2 _ _ /) where rrii counts the occurence of the integer i in 
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the partition (Ai, A2, . . . , \p). Further let ({X) = (HiLi ^"^''^iO ^ and ^(m) denote the 
multinomial coefficient 

- w- («' 

Products of special symmetric polynomials introduced above according to a partition A 
will be written with a subscript A, for example, 

p 

/a^H/^- (20) 

1=1 

Then 

f"" I = E(±1)''^'^'^C(A).A, (21) 
where the +sign refers to the 6's and the —sign to the /'s, and 



= $:(-i)""'^'VM . (22) 



Note the symmetry of the transformations f ^ b and b f. This is a consequence 
of the invariance of the identities ([15D with respect to the interchange b ^ f and of 
^0 = /o = 1- Note that (0), according to ([iTl), is essentially equivalent to a well-known 
explicit expansion formula for the reciprocal value of a power series. These and other 
explicit relations between special symmetric polynomials are often also written in terms 
of determinants, see Ec 



4. Partition functions 

The canonical partition function plays a fundamental role in statistical mechanics since 
most thermodynamic functions can be derived from it. It is defined by 

Z{P) =TTexp{~pH), (23) 

where H denotes the Hamiltonian of the system. Sometimes one writes Zn{(3) in order to 
stress the dependance on the number of particles in the system. The grand canonical 
partition function is defined by 

00 

Z(^,/?)^5^Z^(/3)^^, (24) 

N=0 

where the variable z is physically interpreted as the fugacity of the system. 

Here we will only consider ideal gases, i. e. the case where H is the sum of 
(identical) single-particle Hamiltonians 

N 

H = J2hn, (25) 

n=l 

for example, fermions in a common potential or mean field. 
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Let denote the energy eigenvalues of counted in such a way that several 
Ei have the same value in case of degeneracy. Then the bosonic(+) or fermionic(— ) 
partition functions depend only on the Ee and can be written as 

N 

ZN{f3)= Yl exp(-/?5]i?0 (26) 

il<«2<---<«iV £=i 

or 

N 

ZUf3)= E exp(-/?^i?0. (27) 

ii<i2<...<ijv ^=1 

Let us consider for the moment only systems with a finite number of energy eigenvalues 
El, . . . , En- Expanding the exponentials in (p6D,(p7D into products and introducing the 
abbreviations 

x, = exp{-PEe), e=l,...,n, (28) 
it is easy to see, c. f. section 2, that 

Z^ = fNixu...,Xn), (29) 

and 

Zj^ = b^ixi, . . . ,Xn). (30) 

Hence the fermionic/bosonic partitions functions for particular systems are obtained if 
the corresponding symmetric polynomials fN/i>N are evaluated at the values (^8]) , where 
only (3 is left variable. The grand canonical partition functions obviously correspond to 
the generating functions of special symmetric polynomials introduced in (p!6D . It is then 
obvious that suitable relations between symmetric polynomials can be translated into 
relations between partition functions which are independent of the particular physical 
system, e. g. of the common potential or the dimension of the physical space, as long as 
the systems are ideal gases. 

For systems with infinite-dimensional Hilbert spaces, i. e. n = oo, the polynomials 
degenerate into formal power series. Evaluation at the values (p8D leads to partition 
functions only if these series converge. But then the translation of combinatorical results 
into physical ones is still possible. 

Having identified the physical interpretation of the /'s and the 6's it remains to 
investigate the power sums Sj. Since x\ = exp^—j/SEi) we conclude 

.,(e-^^^...,e-^^^...) = Z±(J/3). (31) 

Note that = Zf . 

Now we can interprete all results mentioned in the last section as relations between 
the physical partition functions. We will not reproduce these relations anew but say 
a few words on each. The Newton identities (|T3p are viewed as recursion relations for 
partition functions which allow the explicit calculation of Z]^(p) if Zi((3) is known for 
all (3. In this context they appear for the first time in the textbook of P. Landsberg 
2T| . Later they have been re-discovered and applied recently to systems of few quantum 
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particles, see e.g. [0, |2y, ^ |^ . An explicit representation equivalent to ( ]21| ) appears 
within a physical context in ^ and Equations (0) and (^) express relations 
between bosonic and fermionic partition functions which have, to our knowledge, not 
yet been considered in the physical literature. Note, however, that the identity ( ]T8| ) can 
be rewritten as 

Z^iz,P)Z-{-z,P) = l, (32) 

which obviously follows from the product representations and is implicitely contained 
in the above-mentioned articles based on a "unified treatment" of Fermi and Bose 
statistical mechanics of ideal gases H, |^. We remark that (|3^) implies = 
which is closely related to the fact that in the classical limit (i. e. 2 — > 0) the distinction 
between Fermi and Bose gases disappears. Further we note that from the product 
representation (0,113) and (|l^) one immediately derives the well-known fact that Z~ 
has zeros at ze = — exp^jSEi) and no poles, and hence has poles at Z£ = exp{PEe) 
and no zeros. Poles and zeros lie outside the physical domain of 2; G (0, 00) for fermions 
and z G (0, 1) for bosons. 

5. Concluding remarks 

In this article we have suggested a new way of looking at the partition functions of 
ideal quantum gases. The construction of these functions is split into two steps. In 
the first step we consider a finite or infinite set of abstract energy levels as symbolic 
variables and define certain symmetric polynomials (or formal power series) in these 
variables: the for bosons and the Jn for fermions. The symmetry involved is a 
symmetry between energy levels, not between particles. It is due to the fact that for 
non-interacting particles all possible occupations of the energy levels are counted on an 
equal footing. This first step is the same for each physical system. 

The second step consists in an evaluation of the symmetric polynomials and /jv 
at the values e~^^^ for the variable X£. Only within this step the numerical energy values 
of the concrete physical systems, their degeneracies et cetera enter into the expression 
for the partition functions. Now the different possible occupations of the energy levels 
contribute to the partition function with different weights. The mentioned physically 
relevant combinatorical results are independent from the second step and hold for all 
considered systems. 

Finally we would like to mention that the wording "partition function" is due to 
Darwin and Fowler and is obviously connected to the "partition" of energy into 



the energy levels of a system. Thus the combinatorical element was present from the 
outset, although the connection to the combinatorical partitions occuring in the explicit 
formulae of section 3 seems to be novel. 
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